A study of elliptic flows in a quark combination model 
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We carry out a detail study of elliptic flows in Au- Au collisions at 200 AGeV in a quark combination 
model. We find that elliptic fiow data for a variety of hadrons can be well reproduced except pions if 
constituent quarks with equal parallel transverse momenta combine into initially produced hadrons. 
In a combination mechanism where initial hadrons are formed by quarks with unequal parallel 
transverse momenta, theoretical predictions agree with data for all available hadrons including 
pions. The mass hierarchy at low transverse momenta in elliptic flows can be understood in the 
same quark combination mechanism as in the mediate range of transverse momenta. 
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I. INTRODUCTION 

The elliptic flow is one of the most important observables in ultra-relativistic heavy ion collisions, which carries a 
lot of information about the initial state of the dense matter Q, Q, |^ 111 . The elliptic flow V2 for a hadron is deflned 
C\j ■ as the coefficient of the term cos 20 in the Fourier expansion of the invariant differential cross section E^^, where (j) 
K>" ' is the azimuthal angle of the hadron momentum in the transverse plane The E877 collaboration flrst measured 

■ the anisotropic flow in heavy ion collisions at AGS energies. In recent years, a thorough experimental study of elliptic 
' flows for various hadrons has been taken at relativistic heavy ion collider (RHIC) energies. The elliptic flows as 
functions of centrality, pseudorapidity and transverse momentum have been measured. Various theoretical models 
have been proposed to describe data, such as coalescence or recombination models 0, @, 0, E|, transport models 
\^ , ^1, 12, 12j, and hydrodynamic models iJJ, 15., .1^] etc.. 

' A salient feature of V2{pt) at mediate transverse momenta pr is the quark number scaling, i.e. when both V2 and 
Pt rescaled by the constituent quark number Uq of the hadron, the eUiptic flow is a universal function of px for all 
hadrons. The quark number scaHng of the eUiptic flow and the abnormal ratio of baryons to mesons at mediate pr 
support the quark recombination mechanism in hadron formation. However, at small transverse momenta pt the 
^ . quark number scaling is not obvious, instead, there is a different structure the so-called mass hierarchy ^3 in data. 
In hydrodynamic models, the mass hierarchy arises from radial flows which are hadron masses dependent. But 
hydrodynamic models cannot describe eUiptic flow data at mediate and high pT- The mass hierarchy structure of 
viip r) for charged pions and kaons and protons / anti-protons has been reproduced by a multi-phase transport model 
[l^ . In the large pt range of about 10 GeV/c, Ref. flil l2(l l2l] | have described the global behavior of V2{pt) for some 
5^ ■ hadrons in combination mechanism, but the mass hierarchy of V2{pt) for various hadrons at low pr has not been 
a ■ fuUy compared with experiment data or reproduced within coalescence/recombination models. Ref. ^,i23] and some 
models have given several possible origins of the mass hierarchy, however, we show decay is one of the most important 
reasons and give another explanation for the mass hierarchy. In this paper, we carry out a detail study of elliptic flows 
in a quark combination model in the whole pr range with a special focus on the mass hierarchy or fine structure at 
low Pt- We will show that the fine structure of elliptic fiows can be well described in the quark combination picture 
provided initially produced hadrons and resonance decays are treated separately and properly. 



II. QUARK COMBINATION MODEL 

All kinds of hadronization models demand themselves, consciously or not, satisfy rapidity or momentum correlation 
for quarks in the neighborhood of phase space. The essence of this correlation is its agreement with the fundamental 
requirement of QCD which uniquely determines the quark combination rule OESl- According to QCD, a qq may 
be in a color octet or a singlet which means a repulsive or an attractive interaction between them. The smaller the 
difference in rapidity for two quarks is, the longer the interaction time is. So there is enough time for a qq to be in a 
color singlet and form a meson. Similarly, a qq can be in a sextet or an anti-triplet. If its nearest neighbor is a i? in 
rapidity, they form a baryon. If the neighbor is a q, because the attraction strength of the singlet is two times that 
of the anti-triplet, qq will win the competition to form a meson and leave a q alone. 
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Figure 1: The V2{pt) spectra for final state anti-proton from Eq. Q (dashed line) and Q (solid line). The sub-diagram in the 
upper- left corner shows curves of Eq. Q (dashed line) and (EJ (solid line). 
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We have developed a variant of the quark combination model (QCM) based on a simple quark combination rule 
[3,12^ incorporating the above QCD requirements. The flavor SU{3) symmetry with strangeness suppression in the 
yields of initially produced hadrons is fulfilled in our QCM |23. E^. Using our QCM, we have described most of 
multiplicity data for hadrons in electron-positron and proton-proton/anti-proton collisions |23. l2^ l2^ 123. lifll 
Also we solved a difficulty facing other QCMs in describing the TASSO data for baryon-antibaryon correlation in 
electron-positron collisions: they can be successfully explained by our QCM (2H . Our QCM can also combine 
with the color fiow picture to describe the hadroniztion of multiparton states I,32L [3^ 0| . We have extended our 
QCM to reproduce the recent RHIC data for hadron multiplicities and px spectra 

In this paper we present a detail analysis of V2[pt) in the whole pr range based on our QCM. We put particular 
emphasis on the low pt region where the elliptic fiow data for various hadrons show non-trivial fine structure. Our 
treatment of elliptic fiows is different from other coalescence or recombination models (lflll2(l |. We extract V2,q{pT) for 
quarks as input from fitting the experiment data of various hadrons. So we can construct elliptic fiows for all initially 
produced hadrons from v^^qipx) for quarks in our QCM. Then we let all resonances decay by using the decay program 
of the event generator PYTHIA j^^l to obtain elliptic fiows for final state hadrons. One advantage of using PYTHIA 
to deal with decays is that all available decay channels are covered. Here we only take into account light quarks in 
hadronization. In central Au-Au collisions at 200 AGeV, the yield of charm hadrons is about 28 in the whole rapidity 
range 37]. So charm quarks are negligible compared to light quarks in V2{pt) for final state hadrons. 



III. RESULTS AND ANALYSIS 



In order to describe data of elliptic fiows, only those quarks with collinear momenta are allowed to combine into 
hadrons, because the combination of constituent quarks in different momentum directions would lead to large invariant 
masses which are off-mass-shell for most hadrons. To describe data in central rapidity, we choose the rapidity range 
y £ [—0.5, 0.5] for hadrons in our calculation. We also assume that quarks and antiquarks of a same fiavor have the 
the same V2.q{pT)- Considering that the coalescence or recombination models work very well at mediate pr, we focus 
our study on elliptic fiows at low pr which have not been well investigated in combination-like models. We assume 
three collinear combination strategies, the naive scheme with equal pr combination and the same V2,q{pT) for quarks 
of all fiavors, the modified naive scheme with equal pr combination and a different V2,q{pT) for strange quarks from 
light quarks, and the scheme with different px combination. 
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Figure 2: The U2(pt) spectra for directly produced (upper panel) and final state (lower panel) hadrons from Eq. 
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A. Naive combination scheme 

In this scheme, all u, d and s quarks are assumed to have the same V2q{pT) spectrum. The quarks with equal pr 
combine into hadrons. Our input is based on 

vf^{pT) = 0.078 tanh(1.59]5T - 0.27) (1) 

for quarks given in Ref. by fitting the proton data. Hereafter the unit of pr is GeV. We make two improvements, 
one is that the decay contributions have been subtracted, the other is that the negative part of v^^ (pr) inconsistent 
to data and to the expected quadratic behavior 113,1131 at small pr (see the dashed line in Fig. 1) is corrected to be 
positive. The improved V2,q{pT) is given by 

V2,q{pT) = 0.256p|, 0<pT< 0.365, 

V2Apt) = 0.078 tanh[1.45(pr - 0.25)], pt > 0.365. (2) 

Note that considering the continuity of the Eq. we adopt the pj< instead of the p^ form. The difference resulted 
from the two function forms at low px range can be neglected for the coefficient of pj^ or p^ is adjustable. One sees 
that the data for antiprotons are well reproduced at small pr with the above input, see the solid Hue in Fig. 
Moreover it can fit all hadron data except pions as shown in the lower panel of Fig. |2l Note that in figures following 
Fig. 12 the order of curves from left to right is the same that of data for various hadrons. The STAR data in Fig. 
are taken from Ref. In Fig. [Hand El the PHENIX data are taken from Ref. In Fig. El the data for K'^ and 

A + A are from Ref. jiJl, the data for S~ + and fl~ + fl^ are from Ref. and the data for kaons are from Ref. 

0- 

In the naive combination scheme, directly produced mesons and baryons have the same elliptic flows V2^m(pt) = 
2i'2,g(pT/2) and V2.b{pt) = 'iv2.,q{PT /i) respectively, regardless of their masses and species. Comparing the upper and 
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Figure 3: The U2(pt) spectra for directly produced (upper panel) and final state (lower panel) hadrons from Eq. 
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lower panels of Fig. |2l one observes that the V2{pt) curve for the hadron i shifts somewhere in the left depending on 
how much the decay contributions the hadron has. One sees that the fine structure of V2(pt) for hadrons at small pr 
results from decay contributions. But for pions it is not enough to reproduce the data by only including the decays, 
one has to take into account other effects. We will address this issue in the modified naive combination scheme in the 
next subsection. 



B. Modified naive combination scheme 



In this scheme, we take into account that V2,q{pT) for strange quarks is different from light quarks. The difference 
is mainly at small pr since the quark mass effects barely infiuence the magnitude of V2(pt) at mediate pt, see the 
attached figure in the upper-left corner in the upper panel of Fig. Note that Ref. [Tflll2(l| also include the strange- 
light difference. We know that the strongly coupled quark-gluon-plasma produced at RHIC is relativistic fiuid, then 
we can assume that the average transverse momenta for u, d and s quarks are proportional to their average energies, 
'^<PT> « Based on this relation, we modify V2,q{pT) for u, d and s quarks as, 

V24pt) = 0.256(pTM)^ 0<ptM< 0.365, 

V2Apt) = 0.078 tanh[1.45(pTM -0.25)], PtM > 0.365. (3) 



where i — u,d,s and at is given by 



2 < ^ >u,d „ r^A 

< ^ >u.d + < -tj >s 

as ^ « 1.12. (4) 

< E >„.<j + <E>s ^ ' 
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Figure 4: The V2{pt) spectra for final state hadrons are given in the upper panel. The «2(pt) spectra scaled by constituent 
quark number are given in the lower panel. The theoretical curves are results in the unequal momentum combination scheme. 
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Figure 5: The «2(pt) spectra for <j> mesons in three combination schemes. 
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where the average values of quark energies can be approximated as {E) ^ « y"^/ + (Pr) with constituent quark 

mass vfif for the flavor / (m„ ^ = 0.34 GeV, = 0.5 GeV) and with (p^) given by the distribution J{pt) = 
{p^^ +Px'^ -\- ° from Ref. Here we have implied the central rapidity y « in estimating the average energy. 
Based on the input for quark elliptic flows in Eq. V2{pt) for meson and baryon are written as 

V2,m{pt) ~ ■i;2,j(PT/2) + W2j(pt/2), 

v2.b{pt) ~ W2,j(pt/3) + W2j(pt/3) + t;2,fc(PT/3), (5) 

with i,i denote the constituent quarks in the meson M and j, k denote those in the baryon B. Therefore the V2{pt) 
curves for hadrons with the same flavor content are identical. The results in this combination scheme are shown in 
Fig. El In contrast to Fig. |21 the V2{pt) curves for directly produced hadrons split according to their flavor contents 
and lead to even larger splittings for final state hadrons. The pion curve is partially improved but not enough to 
match the data. For other hadrons, the agreement between theoretical predictions and data is satisfactory. Although 
tuning au,d to a smaller value would make pion curve fit the data, the fine structure of V2{pt) for other hadrons 
prevent them from splitting too much. The results show that in the quark combination mechanism the quark mass 
effect and resonance decays partially account for the fine structure of V2{pt) for final state hadrons. In contrast to 
the hydro-model interpretation of differences among hadron V2{pt) as hadron mass effects, the mass effects at the 
quark level is partially responsible for such differences in the quark combination picture. To identify the mass effects 
of strange quarks and light quarks, a better way is to measure and compare V2{pt) for (j) and K* in experiments. 

C. Unequal Pt combination scheme 

As shown in Fig. El it is still not enough to reproduce the pion data by including resonance decays and the quark 
mass effect. In order to better fit the pion data together with all other hadrons, we relax in this scheme the constraint 
that quarks with equal pr combine into hadrons by allowing unequal px combination. The unequal pr combination 
is more general than the equal one which can be regarded as a special case of the former. Since the combination 
condition is changed, we then start with a different form of V2.q{pT) for quarks from that in the equal pr combination 
schemes, 

V2,q{PT) ^ 0.286]5:^, 0<pT< 0.525, (6) 
V2,q{pT) = 0.11 tanh[1.6(pT- 0.4)], pT > 0.525. (7) 

The results for hadron V2{pt) are shown in Fig. 4. The initial hadrons with the same fiavor contents have an identical 
curve. One sees in the upper panel that the pion curve agrees with data very well. Note that Ref. j^S] reproduces 
the pion data after taking into account the momentum distribution of hadrons, a similar strategy to ours, and Ref. 
[4^ takes into account the resonance decay contribution to pions. In contrast, we give the results for other hadrons 
also fit the data besides pions. The main difference from the equal pr combination scheme is that V2{pt) spectra for 
hadrons at mediate px split into two groups (of mesons and baryons), i.e. there is no exact quark number scaling, 
but theoretical predictions are still within the error bars. More precise measurements are needed to distinguish the 
two combination schemes in the mediate px range. Obviously the fine structure of V2{pt) at low pr is a joint effect 
of resonance decays, the quark mass difference in V2,q{pT) and unequal pT combination. The reason that the unequal 
PT combination works well at low and mediate pr lies in the fact that it splits the spectra for different hadrons at the 
cost of slightly spoiling the quark number scaling. But the data are not sensitive to such a violation of the scaling. 

In Fig. El we give the predictions for 4> mesons in three combination schemes. The agreement with recent STAR 
data is satisfactory Q|- In all of our calculations, we neglect the infiuence of hadron re-scattering on V2{pt), same as 
in Ref. 

IV. SUMMARY 

We use quark combination model to study eUiptic fiows for directly produced hadrons before resonant decays and 
the final state hadrons. Especially we have investigated the fine structure of elliptic fiows at small pr- We find the 
combination mechanism works very well in describing elliptic fiow data in low pr as well as in mediate pr regime. 
In the scheme with equal pr combination and equal input function for V2,q{pT) for light and strange quarks, the 
data for all hadrons can be well reproduced except pions, where the fine structure can be understood by the decay 
contributions. Taking into account the mass effect of the input function i'2,g(PT) for light and strange quarks and 
assuming unequal pt combination, the overall agreement between our predictions and all available data for hadrons 
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including pions are satisfactory. The decay contributions and the mass effect of the input function U2,q(pr) for quarks 
play important roles. The directly produced hadrons with the same flavor content have an identical V2{pt) spectrum, 
which is an obvious feature of the quark combination mechanism. 
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